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Starting from a multipole expansion of intra- molecular Coulomb interactions, we present con- 
figuration interaction calculations of the molecular energy terms of the hole configurations (/ij) m , 
m = 2 — 5, of CJJ + cations, of the electron configurations t" u , n = 2 — 4, of the Cgg" anions, and of 
the exciton configurations (h^t^ u ), (h^ti g ) of the neutral Ceo molecule. The ground state of Cgg 
is either 3 Ti g or 1 A g , depending on the energy separation between ti g and ti u levels. There are 
three close (~0.03 eV) low lying triplets 3 T lg , 3 G g , 3 T 2g for CgJ, and three quartets 4 T lu , 4 G U , 
4 T2u for CgJ , which can be subjected to the Jahn- Teller effect. The number of low lying nearly 
degenerate states in largest for m = 3 holes. We have calculated the magnetic moments of the hole 
and electron configurations and found that they are independent of molecular orientation in respect 
to an external magnetic field. The coupling of spin and orbital momenta differs from the atomic 
case. We analyze the electronic dipolar transitions (tin) 2 — ► t\ u ti g and (ti u ) 3 — > (ti u ) 2 ii 9 for Cgg" 
and Cgg. Three optical absorption lines ( 3 Ti g — > 3 H U , 3 Ti u , 3 A u ) are found for the ground level of 
Cgo~ and only one line ( 4 A U — » 4 Ti 9 ) for the ground state of Cgg - . We compare our results with the 
experimental data for Cgg" in solutions and with earlier theoretical studies. 

PACS numbers: 31.10.+Z, 31.25.-v, 75.75.+a, 73.21.-b 



INTRODUCTION 



The physics and chemistry of fullerenes is full of sur- 
prises. After almost ten years of intensive theoretical and 
experimental work, unexpected discoveries of ferromag- 
netic polymerized Ceo Q as well as superconductivity 
in hole doped pristine Cgo 0] and in lattice-expanded 
Ceo § raise new questions about the electronic struc- 
ture of the Cgo molecule and its molecular ions. Using 
the field-effect doping techniques, it has been shown that 
high transition temperatures are achieved for the case of 
three holes per the C^o molecule Q , when many electron 
effects come into play. Motivated by these new experi- 
mental findings, here we present a calculation of the elec- 
tronic structure of few holes of C™ + . Our second goal is 
to study the molecular term picture of the Cqq anion 
(n = 2,3,4), which behaves as a quasi-element in a vast 
majority of ionic compounds When C^ ± units are in 
a solid, additional inter-molecular interactions superim- 
posed on the intra-molecular ones should be taken into 
account and a realistic theory of solids should treat both 
kinds of interactions on equal footing. Therefore, the 
consideration of intra- molecular correlations should be 
a necessary ingredient of a many electron theory aiming 
to describe superconductivity, magnetic properties or a 
metal-insulator transition ||, [|, [| [?], ^ . The problem is 
also a challenge from experimental point of view j^, ^) . 
fn the present paper we will give quantitative results on 



the many electron terms and magnetic moments of 
and Cgg". 

While the one-electronic structure of the neutral C 
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molecule is known for many years pdj | (in fact, even be- 
fore the actual discovery of the Buckminsterfullerene | fl2"| ) 
the case of two or three electrons (holes) on the degener- 
ate ti u or h u shell requires a special treatment. Then the 
two (n — 2) or three [n = 3) electrons (holes) are equiv- 
alent and should be treated on equal footing. Such kind 
of situation is very typical for atoms with open electron 
shells. When an atom has two or more valence electrons 
on orbitally degenerate states (like p, d or /) its energy 
spectrum can be very complex reflecting the electronic 
degrees of freedom. There exist empirical observations 
known as Hund's rules, which prescribe the occupation 
of the orbitals, but those are just consequences of the 
atomic theory of many electron states formulated by Con- 
don and Shortley long ago |f3| . The real driving force be- 
hind the term splitting is the Coulomb repulsion of the 
valence electrons. While in atoms electronic energy lev- 
els are split due to the multipolar Coulomb intra-atomic 
interactions, in the case of Cqq we deal with the intra- 
molecular ones. 

Here we present an approach which was inspired by 
the theory of many electron atomic states |l3| . It is also 
an extension of our original method of multipole expan- 
sion for electronic states [ 14| [t5| . It is worth mentioning 
that our treatment should not be confused with an ordi- 
nary single-determinant Hartree-Fock calculation which 
does not take into account the intra-molecular correla- 
tions and the molecular term structure is overlooked. We 
will see (Sec. Ill) that in our approach each basis func- 
tion represents a Slater determinant, and the solution is 
found as their linear combination. As such, the method 
corresponds to a many determinant treatment or con- 
figuration interaction (CI). Therefore, the approach is a 
genuine many electron one as long as we limit ourselves to 
the relevant orbital space (h u , t\ u or ti g ). To our knowl- 
edge, in the literature there are only two calculations 
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concerning the electronic structure of negatively charged 
Cgrj" molecular ions, reported by Negri et al. 



Saito et al. |17 
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and by 

In the latter work, however, the assign- 
ment of excitations with molecular terms is done only for 
the neutral molecule (ra = 0). 

In the present paper we assume that the molecule has 
the icosahedral symmetry. If one wants to consider a 
distortion of or Cg ~" @ H, |j| ||, §§, the com- 
puted energy levels and their eigenvectors can be used 
as a starting point for the description of the Jahn- Teller 
effect in these systems. Indeed, the electron-phonon (or 
vibronic) coupling occurs if [T e i] 2 contains r vi b ^3|, pi[ . 
Here r e ; is the symmetry of the electronic molecular 
term under consideration, while V V ib is the symmetry 
of a vibrational normal mode. It is evident then that 
for a meaningful analysis of the vibrational coupling and 
the resultant Jahn- Teller distortion, one has to know the 
symmetry of the corresponding electronic terms, that is 
the issue of the present study. In addition, now it has 
been realized that the Jahn- Teller effect for CgQ is rather 
weak, and the situation is probably best described as 
dynamic |25|| , where the molecule on average retains its 
icosahedral symmetry. The vibronic coupling for Cg^ + is 
not so well investigated, but estimated as ~ 0.1 — 0.2 eV. 
We will see that the energy span of the hole configura- 
tions a few times exceeds this value. Therefore, if there 
is a static Jahn- Teller effect for C^, + , it can be treated 
in the limit of weak or intermediate coupling, leading to 
a more complex picture where the molecular terms are 
subjected to further splitting 24, 3l| . 

Since the Cgo molecule reveals nontriviai degeneracies 
and peculiarities, its term structure has been studied 
by group-theoretical methods pq, E7L 



29, 30 . The 



analysis there was focused on the symmetry and num- 
ber of terms for g N (N = 1 - 7) f2§ and h N (N=l — 9) 
p?| , |28| , p9[ |30f configurations, but by itself it can not give 
a quantitative picture of the splittings. Where our ap- 
proach overlaps with the group-theoretical one, our find- 
ings are in agreement with the latter. 

The paper comprises the following sections. First (Sec. 
II), we introduce the angular dependencies of h u , t\ u and 
tig functions. Next (Sec. Ill), we describe our method of 
treating multipole Coulomb correlations for many elec- 
trons. Then we apply it to hole- and excitonic configu- 
rations of the Ceo molecule, Sec. IV. In Sec. V we give 
the resultant energy spectra and compare our calcula- 
tions with those of Negri et al. jl6| for Cg^~. In Sec. 
VI we compute magnetic moments, in Sec. VII optical 
lines and line strengths for the electron dipolar transi- 
tions of CgQ and Cg^~. Finally, we give our conclusions 
in Sec. VIII. 



II. ANGULAR DEPENDENCE OF 
MOLECULAR ORBITALS 

The neutral Ceo molecule has the highest molecular 
point group (Ih)- Expanding its density in multipole se- 



ries one finds that nontriviai angular dependencies are 
given by the symmetry adapted functions (SAFs) |3^] 
of A\ g symmetry with the lowest components character- 
ized by I = 6, I = 10 and I = 12 ||, ||. That makes 
the Ceo fullerene the most spherical molecule among the 
others. Owing to such unique symmetry, the parent- 
age of 7r-molecular orbitals in spherical harmonics can 
be clearly traced jy]]. For lowest occupied 7r-levels of 
the neutral molecule we have a g (I — 0), t\ u (l = 1), 



h g (I = 2), t 2u (I = 3), g u (I = 3), g g (I 



4), h g (I 
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and h u (I = 5). The corresponding molecular states ac- 
commodate 60 electrons and the seven electron shells are 
completely filled. Then the generalized Unsold theorem 
Q ensures that the resulting electron density of the tt- 
states has full (or unit) icosahedral symmetry A\ g . (Here 
and below we use capital letters for the irreducible rep- 
resentations (irreps) of density and small letters for the 
irreps of electron wave functions.) 

Due to the direct correspondence between these 7r- 
shells and the molecular orbital index I, one can imme- 
diately find out an orbital part for a given 7r-state in the 
same way as we know the orbital parts of s, p, d and 
/ electrons in an atom. Therefore, the type of an irre- 
ducible representation (a, t\, t 2 , g and h) and the orbital 
index I uniquely determine the angular dependence of 
the molecular orbitals. The angular functions are called 
spherical harmonics adapted for the icosahedral group 
Ih- Such symmetry adapted functions (SAFs) were tab- 
ulated by Cohan in Ref. [52] for all I < 14. (Unfortunately, 
Cohan worked with unnormalized spherical harmonics so 
that it requires some efforts to express SAFs in conven- 
tional spherical harmonics Y™.) In particular for the 
highest occupied molecular orbital (HOMO) h u (I = 5) 
one has 



ipi{h u ) = Yg 5 ' 5 , 




(2.1a) 
(2.1b) 

(2.1c) 

(2.1d) 

(2.1e) 



Here the normalized real spherical harmonics are defined 
with the phase convention of Ref. 31 and the orientation 
of the Ceo molecule corresponds to the choice of the z 
axis as one of 12 fivefold axes and the y axis as one of the 
twofold axes perpendicular to z [^2| . We call this position 
of Ceo the orientation of Cohan. In order to transform the 
molecule to the standard orientation |35]| where molecular 
twofold axes lie along the Cartesian x, y and z direction 
one has to rotate the molecule anti clockwise about the 
y axis by an angle f3 = arccos (2/V10 + 2V5) w 58.28° 
]l4j . Each of the h u orbital function then is expressed in 
terms of Y§ , where r stands for m 



or (m, c), (m, s) 
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of real spherical harmonics (Appendix B). In the follow- 
ing we will work with the molecule in the orientation of 
Cohan. (Of course, the results are independent of the 
choice of the molecular orientation.) 

The lowest unoccupied molecular t\ u orbital (LUMO) 
corresponds to Z = 5 and has the following three angular 
components in Cohan's orientation of C§$: 



(2.2a) 
(2.2b) 
(2.2c) 



The angular parts of t\ u LUMO have been derived before 
in Refs. p36] and 14, In the standard orientation of Ceo 




they are given by Table I of Ref. The LUMO-HOMO 
energy gap is about 2.7 eV §, |,]37j. 

Finally, at an energy about 1.15 eV [38| above 
LUMO one finds the molecular t\ g level with I = 6 
(LUMO+1). In the orientation of Cohan the angular 
parts are given by 



Mhg) = Y 6 
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(2.3a) 
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(2.3b) 
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(2.3c) 



III. METHOD OF CALCULATION 

Our method of multipole expansion of the Coulomb 
interaction has been reported before |lq| . Here we 
extend it and apply to the case of the Ceo molecule with 
the icosahedral symmetry. In the following we consider 
in detail the case of two and three t\ u electrons. Starting 
with a pair of electrons we will give a special attention 
to the procedure of adding one extra t lu electron to the 
pair. In the same way one can add a fourth electron to 
the group of three electrons and etc. Therefore, our main 
goal of treating n electrons can be reached by adding one 
electron after another. 

Since the estimated one-electron spin-orbit coupling 
is negligible (~ 0.16 cm -1 ) fl36| , we are working in the 
a LS (Russell-Saunders)" molecular approximation. (The 
spin-orbit coupling is a single particle operator and in 
principle can be included in the calculation [jl5).) We 
start with a pair of t± u electrons and label the two- 
electron basis ket-vectors by a single index I which in- 
corporates a pair of one-electron indices (21,22)1 



\h; 12) 



(3.1) 



The indices i = (k,s z ) stand for the t\ u orbitals (k = 
1, 2, 3) and the spin projection quantum number. The 
corresponding basis wave functions are 



(r,r>\I) = {r\h) ■{?%), 



(3.2) 



where (r \i) = lZ(r)(n\i) . Here 1Z is the radial component 
of the t\ u molecular orbitals (MO), n stands for polar 
angles Q = (O, 4>). There are six orientational t\ u vectors 
(or spin-orbitals) (n\i) (i.e. i = 1 — 6), 



(n\i) = ip k {n) U s (s z ). 



(3.3) 



Here ipk are the three t\ u MOs as given by Eqs. ( [2.2a -c) 
for the Cohan's orientation of Ceo , u s is the spin function 
(s = ±) for the spin projec tion s s z = ±1 /2 on the z-axis. 

The order of indices in ( |3.l| ) and (3.2) is important if 
we associate the first electron with the state i-y while the 
second with the state 12 ■ From the dynamical equivalence 
of the electrons we can permute the spi n-or bitals of the 
state \%2] ii) to the standard order, Eq. (|3.lD, by using 



*2; n) 



In; 12) 



(3.4) 



since it requires the interchange of the two electrons. In 
order to descr ibe the same quantum state we w ill u se the 
basis vectors (3.1) where i\ > i 2 and apply Eq. ( [L4| ) when 
needed. (Alternatively, one can use the standard proce- 
dure of antisym mctr ization of the basis vectors (3.1).) 
Thus, our basis (3J) consists of (6 x 5)/2 =15 different 
vectors \I). 

In the following we will study the intra- molecular cor- 
relations of electrons within a formalism based on a mul- 
tipole expansion of the Coulomb potential between two 
electrons (charge e = — 1), 



V(r,f') = 



1 



(3.5) 



The multipole expansion in terms of real spherical har- 
monics Y®, Y™ ,c and Y™' s (we use the phase convention 
and the definitions of Ref. [Bll) reads: 



V(r,r') = J2Mr,r')Y l T (n)Y l T (ti), 



(3.6) 



where t stands for m = 0, (m, c) or (m, s) of the real 
spherical harmonics and 



vi(r, r') 



4tt 



„(*+!) ) 21 + 1' 



(3.7) 



with r> = max(r,r'), r < = min(r,r'). 

The direct matrix elements for the intra-molecular 
Coulomb interactions are obtained if we consider the 
i\ ~ * ji transitions for the first electron and the 12 — > 32 
transitions for the second (we recall that i\ > i 2 and 
ji > 32). We label this two-ele ctro n transition by the 
index 02 = 1. Starting from Eq. (3.6) we obtain 



</|U(r,r')|J) Lou ' =^v l c l , r {i l j 1 )c l , T {i2j2), (3. 

l,r 
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where 



vi = Jdrr> /drV'fcVJttVMry) ( 3 - 9 ) 

accounts for the average radial dependence. The transi- 
tion matrix elements ci tT are defined by 

c ltT (ij) =JdQ (i\n) Y{{h) (n\j). (3.10) 

The other possibility is to consider the transitions i\ — > 
j 2 for the first electron and the transitions i 2 — ► ji for 
the second. We label it by the index a 2 — 2. This gives 



the exchange interaction and then we use (3.4) to return 
to the standard order of the spin-orbitals. We find 

(I\V(f,f')\J) exch = -^1^132)0^231), (3.11) 

l,T 

where vi again is given by Eq. ( [3.9] ) and the coefficients 
ci :T by Eq. ( 3.10 ). We observe that in the basis with the 
real t\ u orbitals, and with the real spherical harmonics 
Y{ the coefficients c/ iT are real. 

We start with the spherically symmetric term (I = 0) 
corresponding to th e triv ial function Y° = 1/^/4tt. The 
coefficients q iT in (3.10) become diagonal, c; = o(ij) = 



l/V47T(5y. In considering the other contributions (with 
I > 0) we take advantage o f the selection rules imposed 
by the coefficients c/ iT , Eq. ( 3.1C| ). First of all, we notice 
that the coefficients c; r arc diagonal in terms of spin com- 
ponents u s . Secondly, the odd values of I are excluded 
due to the parity of the integrand in ( 3.10| ), and from 
the theory of addition of angular momenta we know that 



I, 



10. From these observations we conclude that 



I = 0, 2, 4, 6, 8 and 10. Earlier the coefficients c; jT have 
been used for the descri ption of the crystal field of C^~ 
in Ref. @. By adding ( |3.1l| ) to (|J) we obtain 



(I\V(r, r')\J) = U 6(1, J) + v 2 c 2 {I\J) + v iCi {I\J) 
+v 6 c 6 {I\J) + v 8 c s {I\J) + v w c w (I\J), (3.12) 

where Uq = vq/Att is the Hubbard repulsion, 6 the Kro- 
necker symbol, and 



ci(I\J) 



E 



C(Ur){hjl) C(Lr){i232) 
C(Lr){hj2) C{l iT ){l2h)]- 



(3.13) 



We have studied the secular problem for the 15 x 15 ma- 
trix of intra-molecular interactions and obtained 15 en- 
ergy levels E p [t\ u ] (p = 1 - 15), 



Ep[*l«] = Uq + v 2 X 2 {p) + v 4 X 4 (p) 
-H> 6 A 6 (p) + v s X s (p) + wioAio(p), 



(3.14) 



where A; (p) are numerical constants (called "integral" or 



"molecular invariants" in Refs. 126 



(ca 

s 



are quoted in Table I. The 15 levels of (3.14) form three 
distinct terms, i.e. a 9-fold degenerate {t{ u ; 1}, a 5-fold 



They 



TABLE I: Coefficients \[{p) = \i{p) x 10 3 for (t lu ) 2 (p = 



15) and (ti„) 3 (p' 



20). 





p 


deg. 


/ = 2 


4 


6 


8 


10 




1-9 


(9) 


-2.943 


-21.190 


21.817 


-39.706 


-37.556 




10-14 


(5) 


0.589 


4.238 


23.596 


7.941 


11.641 




15 


(1) 


5.886 


42.379 


26.265 


79.413 


85.436 




1-4 


(4) 


-8.829 


-63.568 


65.450 


-119.119 - 


112.667 


K 


5-14 


(10) 


-3.532 


-25.427 


68.119 


-47.648 


-38.872 




15-20 


(6) 


0.0 


0.0 


69.898 


0.0 


10.324 




Q 3 =l, P 


3= + l 




« 3 =2, P 3 


=-1 


a 3 =3, P 3 


=+1 




classification ,„ „ v 
(2,1) index (2,2) 

P=-l parity P=+l 



FIG. 1: Diagram for the calculation of a matrix element 
{111213 ) by expressing it in terms of different one- 

electron transitions i — > j. See text for details. 



degenerate {t\ u ]2} and a single-level term {tf M ;3}. The 
symmetry between the four electrons and the two holes 
within the t\ u LUMO configuration implies Xi(tf u ; p) = 
Xi(t\ u - p), p = 1 - 15. Therefore, Eq. (pl^) holds also 



for the case of four t\ u electrons if we write 6t/o in place 
of f/o- In order to study the splitting quantitatively, one 
has to calculate the radial integrals vi, Eq. (p\9|). We 
leave this calculation for next section and now move on 
to the case of three t\ u electrons. 

For the t\ u there are 6x5x4/3! = 20 basis ket- vectors 
|/') = 1*1*2*3) 5 where i\ > 12 > «3- Now for the Coulomb 
interaction we have a sum of three two-body terms, 

=V(r 1 ,r 2 ) + V(r 1 ,r 3 )+V(r 2 ,r 3 ), (3.15) 

where each V(f a , ff,) is given by the multipole expansion 
( |3 . 6[ ) . In considering a matrix element (/'|I<^ 3 ) |J') we 
have many subcases which we also call transitions. We 
sort them out as shown in Fig. 1. It is important to 
notice that for the third electron there are only three 
possibilities, i.e. i 3 — > j'3, 13 — > j 2 and 13 — > j\, which 
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are labeled by a 3 = 1,2 and 3, respectively. From the 
anticommutation relations we find that the parity (P3) 
of the transitions are +1, —1 and +1. Each case leaves 
possibilities for further transitions of the remaining two 
t\ u electrons. We denote their final states as j[ and j' 2 . 
(For 03 = 1 these indices are j± and j'2; for 03 = 2 j\ 
and j3, Fig. 1; and for 03 = 3 ji and j'3.) However, 
the two-electron transitions have been sorted out before 
when we considered the t 2 u configuration. In that case 
there are only two options, (02 = 1) ii — > j[, «2 - * i 2 i 
and {0,2 = 2) £1 — > j 2 , *2 - * Ji- The first has parity 
P 2 = +1, the second P2 = — 1. 

Therefore, we can classify the transitions by two in- 
dices (03,02) where 03 = 1 — 3 and a 2 — 1,2. For a 
general case of n electrons it is [a n , a>n-i, d n -2, 0,2), 
where = 1,2, ...,k. The parity of an (03,02) tran- 
sition is P = P 3 {a 3 ) P 2 (o2), for the general case P = 
Pn(^n)Pn— i(ttn— 1) ' ••• ' p2(o2). The total number of tran- 
sitions of three t lu electrons is 3 x 2 = 6, in general it 
is nl. For each of the subcases we calculate the matrix 



element of V^ 3 \ Eq. (pjl) 



{l'\V^W a ^ = P 3 (a 3 )P 2 (a 2 ) 

X ^ Q tT (iij'i ) C;, T (t 2 Ja ) ^fe ) ja 3 ) + P-*- ( 3 -!6) 



Here p.i. stands for the other pair Coulomb interactions 
V(r a , 7%). For three particles there are thre e terms, Eq. 
( 3.15 ). The first term is given in Eq. (|3.16|) , two others 
arc found by replacing cj )T (iij[) ci. T {i 2 j 2 ) 8{fa,3a 3 ) witn 
ci, T (iij[) ci yT (i 3 ja 3 ) S(i 2 ,j' 2 ) and with ci >T (i 2 j 2 ) ci lT (i 3 j aa ) 
8(ii, (The parity of the transitions of course remains 
the same.) For the general case of n electrons one has 
(2) different pairs for each n-electron transition. (Here 
(2 ) is a binomial coefficient.) Finally, the matrix element 
\r\VW\J') is found as 



(7'|y(3)| j') = £ Y,{I'W^\J') 



(0.3,0,2) 



(3.17) 



03 = 1 a 2 = l 



Having found the matrix elements, we solve a 20 x 20 
secular problem for t\ u . The energy levels are given by 
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EAt 3 lu ] = ( n f)U + \vMp') 



1=2 



(3.18) 



where only even values of / occur, p' = 1 — 20, and A; 
are quoted in Table I. There are three terms, a 4-fold 
degenerate {<?„; 1} (ground state), a 10-fold degenerate 
{£?„; 2} and a 6-fold degenerate {t? u ; 3}, see Table I. 

The classification scheme described above is very useful 
for handling a single-particle interaction A. In particu- 
lar, the electron coupling to an external magnetic field, 
the spin-orbit interaction, crystal electric field effects and 
the electron dipolar operator fall in this class. The main 
difference from the Coulomb case is that now the interac- 
tion occurs to a single electron (represented by an arrow 
in Fig. 1) while the rest (the other (n — 1) electrons, or 



arrows in Fig. 1) produce (n — 1) Kronecker factors. For 
example, for the three t\ u electrons we have 



A = Ax + A 2 + A 3 , 



(3.19) 



where A^ refers to a single electron. For each three- 
electron subcase (03,02), Fig. 1, we obtain 

(J'|A|J') (03 ' 02) =P3(03)P2(a 2 ) 

xY,(n\A 1 \j[)5(t 2 j' 2 )5(i 3 ,j a3 )+ c.p.. (3.20) 

l,T 

c.p. stands for the two other terms obtained from the 

first by two cyclic permutations, i.e. when zi — > Z2 — ^ i 3 

and j[ — > j 2 — > j a3 . For the general n-electron case 

there will be n such terms for each matrix element 
(7'mj')(«».o«-i.-,a 2 )_ 

Finally, we would like to mention that an operator R 
(rotation or inversion) of the icosahedral group Ih acts 
on all electrons simultaneously and therefore can be writ- 
ten as R = R n R n -i ■ ... ■ R%, where Rk stands for the 
corresponding one-particle operator. The classification 
scheme again is useful for determining the matrix ele- 
ments of the transformation in the many electron space. 
For example, for three t± u electrons we obtain 



3,02) 



(3.21) 



03 = 1 a 2 = l 



where 



mj) (a 3 ,a 2 ) = P^^R^^R^^R^). 

Expressions of that type were used to calculate the char- 
acters of molecular terms and to identify their symme- 
try. The orbital part of the many electron wave func- 
tion transforms as an irrep of Ih , while the spin function 
as a single {t\ u ) or double valued it\ u ) representation 
of SO(2>). Therefore, we classify the molecular terms 



26 21 28 



29| by the symbol 25+1 T, where 25 + 1 is 
the spin multiplicity and T is an irrep of Ih- Thus, the 
molecular terms are 3 Pio, 1 H g , 1 A g for t\ u , and 4 A U , 
2 H U , 2 T lu for t?„, Table | 



IV. APPLICATION TO OTHER 
CONFIGURATIONS 



Below we apply our method for the hole configurations 
(/i+) m of C^ + and excitonic configurations h+t^ u , h^t^ g 
of the neutral molecule. We also consider configurations 
(tiutig) and (ti u ) 2 tig, because they are important for cal- 
culations of electronic dipolar transitions {t\ u ) 2 — > t\ u t\ g 
and (ti u ) 3 — > {t\ u ) 2 ti g of the anions G 2 ^ and Cg^~. 



A. Hole configurations (h^ 



(m - 



5) 



The treatment of a m-hole configuration (h^) m for- 
mally coincides [Ol with the analogous electronic case 



6 



TABLE II: Coefficients m(p) x 10 3 for (hf) 2 . The coefficients 
marked by (*) are not unique depending on integrals vi, see 
text for details. 





deg. 


I = 2 


4 


6 


8 


10 




(9) 


-25.702 


-4.238 


-8.063 


-30.424 


-11.151 


3 G g 


(12) 


-18.835 


-4.238 


-35.345 


-2.015 


-19.244 


Z T 2g 


(9) 


-10.987 


-39.553 


11.212 


-28.746 


-11.503 


Gg 


(4) 


-18.050 


0.471 


29.597 


13.426 


11.721 


1 4 


(1) 


74.162 


59.330 


57.274 


74.228 


136.987 


* 'Hg 


(5) 


-18.492 


0.471 


-12.595 


29.970 


6.970 


* ^g 


(5) 


55.180 


16.951 


8.429 


-18.442 


64.255 



TABLE IV: Coefficients m(p) x 10 3 for {hit^J. 





deg. 


I = 2 


4 


6 


8 


10 


'"Hg 


(15) 


-6.474 


21.189 


-0.381 


11.652 


-23.332 


\H g 


(5) 


42.967 


60.743 


17.919 


16.694 


18.179 


Gg 


(12) 


-8.240 


0.0 


11.184 


-10.435 


19.526 


Gg 


(4) 


-8.240 


14.126 


13.725 


24.304 


103.205 


S T lg 


(9) 


6.474 


-21.189 


9.701 


-11.652 


38.081 


'Tig 


(3) 


6.474 


-21.189 


93.578 


-11.652 


70.528 


3 T2g 


(9) 


15.303 


-14.126 


-0.678 


6.145 


11.644 


'Tig 


(3) 


15.303 


-14.126 


-0.678 


70.660 


69.707 



TABLE III: Coefficients m(j>) x 10 3 for {Jiff. The coeffi- 
cients marked by (*) are not unique depending on integrals 
Vi, see text for details. 







deg. 


I = 2 


4 


6 


8 


10 




4 G u 


(16) 


-55.916 


-33.903 


-61.552 


-39.129 


-48.232 




4m 
-t lu 


(12) 


-62.783 


-33.903 


-34.370 


-67.538 


-40.139 




4rrn 

J-2u 


(12) 


-48.068 


-69.218 


-15.095 


-65.860 


-40.491 




2 T lu + 2 


T 2u (12) 


-0.392 


-5.650 


-30.218 


-13.676 


30.095 




2 - / 
J2« 


(6) 


-62.194 


-5.650 


18.710 


-11.293 


-3.192 




2rp 

-t lu 


(6) 


-47.480 


-40.966 


37.985 


-9.615 


-3.544 


* 


Z H U 


(10) 


-60.006 


-41.550 


-43.897 


-50.436 


-12.533 


* 


2 H U 


(10) 


-51.236 


-5.650 


-24.355 


-5.511 


-10.208 


* 


H u 


(10) 


28.220 


6.121 


-1.637 


24.146 


61.563 


* 


2 Hu 


(10) 


45.941 


44.376 


58.833 


54.239 


108.509 


* 


G u 


(8) 


-2.237 


-40.966 


-36.070 


-43.844 


38.009 


* 


G u 


(8) 


38.730 


6.121 


10.229 


-40.794 


60.323 



(h u ) m , i-e. the Coulomb repulsion between holes and be- 
tween electrons is the same. The only difference concerns 
the spherically symmetric term (Hubbard Uq). Below we 
count the energy from the level of an empty h u shell. In 
constructing the basis functions, one should take into ac- 
count that the one-hole index if changes from one to 
ten, where the five- fold degeneracy is due to the orbital 
freedom, Eqs. ( 2.1a -e), and the two- fold degeneracy due 
to its spin. As a result we obtain the term energies 



+m] /10— m\ 



E P [h+ m ] = ( 



Uo 



10 

Vljl] 

1=2 



(P), 



(4.1) 



where only even values of I occur. The coefficients 
fi[ (p) are again molecular invariants. We quote them in 
Tables [0] and III for two and three holes, configurations 
(h+) 2 and (h^y. Some coefficients depend on vi. This 
did not occur to t\ u or t\ ul Table I. These coefficients 
are marked by star (*) in Tables [n] and III. For them 
(//*) we give values which are calculated with only one 
parameter vi. For example, [i\ corresponds to the case 
when vi = 2 ^ 0, while the others are zero, vi^2 — 0, and 
etc. Interestingly, the stared terms have off-diagonal ma- 
trix elements in the approach of Plakhutin et al. (Tables 
2 and 3 of Ref. 30). The appearance of /i* implies that 
the computation of energy splittings, Eq. (|4.1[), can not 



be separated in two independent evaluations of vi and 
Hi. Finally, we remark that the "accidental" degener- 
acy of 2 T\ U and 2 T2 U states of (h u ) 3 has been thoroughly 
studied in Refs. 



27 



29, B0 



B. Excitonic configurations hut lu and hft 



In order to describe the excitonic configuration (hft^ u ) 
we introduce one h u hole (missing electron in HOMO) 
and one t\ u electron. The basis functions read 



\I) 



l hu> l uh 



(4.2) 



where now there are 10 states of the h u hole (index z£ ) 
and 6 states of the t\ u electron (index The total 

number of basis functions is 60. The important thing 
here is that we have to treat exchange differently. If we 
consider I as the initial state and | J) = \j^ u ; j^) as a fi- 
nal state, then the exchange transition I — > J is described 
as two electronic transitions %u — * ihu and j u — ► jhu- 

For 

the direct Coulomb interaction we consider i u — > j u and 
jhu —> ihu- In addition, the sign of the direct Coulomb 
and exchange interactions has to be reversed [jl]|. That 
is, 



(4.3) 



for the direct Coulomb (even I) and exchange matrices 
(even I for (hfti u ) and odd I for (h^t^ )). The resulting 
energy spectrum for both configurations is given by 



E p = Ae - 



1 



(4.4) 



where p = 1 — 60. Here the energy of the closed shell 
(h u ) 10 is taken as zero and Ae is an energy associated 
with the promotion of one electron from the h u level to 
the t lu (Ae =2.7 eV §, |||) or t Xg shell (Ae =3.85 eV 
U Q) (see also Sec. V). The calculated values of vi(p) 
are quoted in Table IV for (hf t t^ u ) and in Table |v| for 
(hfti g ). The coefficients vi(p) in the Tables correspond 
to u, with the plus sign, i.e. w ; = vi > 0, where vi are 



given by Eq. (3.9) 
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TABLE V: Coefficients i/j(p) x 10 3 for (h+t 







I = 2 


4 


Q 


t s 


10 




3 G„ 


(12) 


21.425 


0.0 


-19.130 


-5.219 


18.167 




1 G„ 


(4) 


21.425 


0.0 


-19.130 


-5.219 


18.167 




1 1 u 




16.834 


-16.204 


20.434 


4.128 


12.320 




11 U 




16.834 


-16.204 


20.434 


4.128 


12.320 




J- lu 




-16.834 


16.204 


-14.058 


-4.128 


-27.802 




J- lu 


CD 


-16.834 


16.204 


-14.058 


-4.128 


-27.802 




3 To 

lu 




-39.789 


10.802 


21.449 


4.207 


-55 659 




J- 2u 


CD 


-39.789 


10.802 


21.449 


4.207 


-55 659 








I = 1 


3 


5 


7 


g 


11 






0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


G u 


(4) 


0.0 


0.0 


0.0 


19.382 


84.172 


25.116 




(15) 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


H„ 


(5) 


0.0 


0.0 


0.0 


45.811 


3.826 


34.493 


3/T"i 
i lu 


(9) 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1 Tiu 


(3) 


91.820 


0.0 


83.877 


13.282 


0.261 


63.149 


3/T"i 
J-2u 


(9) 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1 T2u 


(3) 


0.0 


142.831 


21.746 


5.520 


0.662 


106.374 



TABLE VI: Coefficients A;(p) x 10 3 for (t lu t lg ); A 3 = 0. 





deg. 


I = 2 


4 


6 


8 


10 


l A 


(1) 


-15.303 


-32.407 


5.797 


-22.605 


-15.059 


3 A 

s±u 


(3) 


-15.303 


-32.407 


5.797 


-22.605 


-15.059 


3 H U 


(15) 


-1.530 


-3.241 


14.927 


-2.260 


-4.757 


1 H u 


(5) 


-1.530 


-3.241 


14.927 


-2.260 


-4.757 


J lu 


(9) 


7.652 


16.204 


21.014 


11.302 


2.111 




(3) 


7.652 


16.204 


21.014 


11.302 


2.111 




deg. 


I = 1 


5 


7 


9 


11 


L A 


(1) 


-64.274 


41.939 


-13.747 


31.956 


4.127 


3 A 

-flu 


(3) 


64.274 


-41.939 


13.747 


-31.956 


-4.127 




(15) 


-32.137 


-4.194 


-8.442 


-14.152 


-23.791 




(5) 


32.137 


4.194 


8.442 


14.152 


23.791 


1 lu 


(9) 


-32.137 


-20.969 


-9.487 


-34.238 


-41.028 


Tin 


(3) 


32.137 


20.969 


9.487 


34.238 


41.028 



C. (tiutig) configuration 



The basis functions here are 



(4.5) 



where, as before, the index i u stands for the t\ u LUMO- 
level (i u — 1 — 6), while the index i g — ( k, s z ) stands 
for the three t lg (LUMO+l)-orbitals, Eqs. fl2.3a -c), and 
the spin projection s z . Thus, i g = 1 — 6 and in total, 
t her e are 6 x 6 = 36 nonequivalent basis vectors |/), Eq. 
(3.1). From the theory of addition of angular momenta 
and selection rules we deduce that for the direct Coulomb 
interactions only the even values I = 0, 2, 4, 6, 8 and 10 
are relevant, while for exchange these are odd numbers 
from one to eleven. The calculated A;(p) are quoted in 
Table 0. 



TABLE VII: Coefficients X t (p) x 10 3 for {t lu )H lg ; A 3 = 0. 
The coefficients marked by (*) are not unique depending on 
integrals vi, see text for details. 



deg. 



/ 



6 



10 



*H g (20) 
2 G a + 2 T 2a (14) 
(12) 
(4) 
(2) 



is 



2 A B 



-1.413 
-2.472 
-10.595 
12.360 
12.360 



-17.949 
-2.244 
-37.393 
11.218 
11.218 



54.714 
53.450 
48.627 
63.844 
63.844 



-37.446 
3.420 
■51.009 
-17.102 
-17.102 



-43.636 
2.127 
-50.504 
-33.333 
-33.333 



2 Ti g 
2 T lg 
2 Ti g 



(10) 
(10) 

(6) 
(6) 
(6) 



-1.413 
11.301 
-17.976 
13.738 
-1.413 



-17.949 
26.923 
-30.048 
53.979 
-17.949 



54.714 
62.580 
42.958 
63.568 
54.714 



-37.446 
23.765 
-12.337 
83.868 
-37.446 



-43.636 
12.429 
-4.859 
79.472 
-43.636 



deg. 



I = 1 



7 



9 



11 



*H g (20) 
2 G g + 2 T 2g (14) 
4 Ti 9 (12) 
(4) 
(2) 



-64.274 
-32.137 
0.0 
-64.274 
32.137 



-16.775 
-4.194 
■41.939 
-41.939 
20.969 



-17.407 
-8.442 
-2.614 

-18.975 
9.487 



-38.347 
-14.152 
-50.216 
■68.476 
34.238 



-56.201 
-23.791 
-43.092 
-82.056 
41.028 



2 H g 
2 Ti g 
2 Ti g 
2 Ti g 



(10) 
(10) 
(6) 
(6) 
(6) 



-32.137 
32.137 
-64.274 
-32.137 
96.411 



-21.987 
13.599 
-41.939 
-17.803 
34.578 



-9.260 
8.738 
■14.929 
-9.006 
22.366 



-33.523 
23.480 
-33.480 
-26.724 
30.075 



-39.202 
30.584 
-40.641 
-18.081 
32.867 



D. Three electron configuration Cti u ) 2 ti g 

In case of the three-electron (ti M ) 2 ti g configuration, we 
construct (6 x 5/2) x 6 = 90 basis vectors 



I/') 



(4.6) 



Here i u \ and i u i are indices of the t\ u LUMO states, i.e. 
iui,iu2 = 1 — 6 and i g labels six t\ g states (i g = 1 — 6). 
Since in that case we deal with two equivalent t\ u elec- 
trons, the basis functions are taken with i u \ > i u i . The 
calculated A; (p = 1 - 90) are quoted in Table |Vl|. The 
important peculiarity of A/ is that as in the case of few 
holes, Sec. IV-A, the coefficients marked by star (*) in 
Table VII (for two 2 H g and three 2 T U terms), exhibit de- 
pendence on vi. For these coefficients (A*) we give values 
which are calculated with only one parameter v\ . For ex- 
ample, Aj corresponds to a calculation where vi = 2 =/= 0, 
vi^2 = and etc. 

Another very interesting observation is that the energy 
accidentally" coincides with the 2 T„ 



G g term 



of the 

states, Table [VII 
the ti u (ti g ) 2 configuration 



The same feature has been found for 



V. ENERGY LEVELS 

In order to study the splitting quantitatively we cal- 
culated the integrals vi using three models for radial de- 
pendence TZ of tiu MOs. In the first model (I) we assume 
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DC 



1 1 1 1 1 1 1 


i 1 i 1 i 1 

\ 


^ *" "i ' i 1 i 1 


1 ■ 1 ■ 1 1 



1 2 3 4 5 6 7 

radius 



TABLE IX: Molecular terms and their degeneracies (in 
parentheses) for (ti u ) 2 and (ti u ) 3 calculated with models I, II 
and III; in eV. (2)^0 is zero of energy. 





3 T lg (9) 


'H a (5) 


l A B (1) 


4 A U (4) 


(ti u ) 3 
2 H U (10) 


2 T lu (6) 


I 


-0.275 


0.175 


0.848 


-0.824 


-0.150 


0.299 


II 


-0.122 


0.077 


0.375 


-0.366 


-0.068 


0.131 


III 


-0.117 


0.073 


0.359 


-0.351 


-0.066 


0.125 



FIG. 2: \TZ(r)\ for model II (dashed line) and model III (full 
line) as a function of r (in A). 



TABLE VIII: Calculated v t for models I, II and III; in eV. 



model 


I = 2 


4 


6 


8 


10 


12 


I 


10.195 


5.664 


3.921 


2.998 


2.427 


2.039 


II 


6.919 


3.064 


1.752 


1.137 


0.797 


0.589 


III 


6.798 


2.965 


1.673 


1.074 


0.747 


0.550 




I = 1 


3 


5 


7 


9 


11 


I 


16.991 


7.282 


4.634 


3.398 


2.683 


2.216 


II 


13.145 


4.390 


2.270 


1.395 


0.944 


0.681 


III 


13.018 


4.279 


2.181 


1.324 


0.889 


0.637 



a delta dependence, i.e. 

K(r)=5(r-r Cm )/r 2 , 



(5.1) 



where rc 60 — 3.55 A is the radius of the Cgo molecule. 
This gives 



v, 



-in 



21 + 1 r Ceo 
In the second model (II), Fig. 2, we use 

Tl(r) = C*exp(-V2|^| 



rc 6 



I), 



(5.2) 



(5.3) 



where C is determined from the normalization condition 
and E — —5.863 eV, the energy of the carbon p z atomic 
orbital in atomic calculations in local density approxima- 
tion (LDA) (such dependence corresponds to the large 
distance limit for the carbon j^-orbital) . In the third 
model (III) we take 



K(r) = C'KpMr - r Ce 



(5.4) 



where C is a normalization constant and TZ Pz is the ra- 
dial dependence of the carbon p z orbital in LD A, Fig . 2. 
The calculated value s of v i are q uoted in Table VIII. 

First, from Eqs. ( B.lj ) and flgjj ) we calculate the 
molecular terms of i{ u and t\ u , Table [X. We observe 



that despite their differences, models II and III give very 
close values. Therefore, we believe that the calculated 
parameters vi of models II and III are realistic and will 
change little if a more refined calculation of 1Z is made. 
On the other hand, the first model is a rude approxi- 
mation, and in the following we will not use it. Notice 



TABLE X: Molecular terms of {t lg ) 2 and (ti 9 ) 3 calculated 
with models II and III; in eV. (2 )Uo is zero of energy. 





3 Ti 9 
(9) 


(U 9 f 
X H a 
(5) 


X A g 
(1) 


4 A a 
(4) 


(ii 9 ) 3 
2 H g 
(10) 


2 Ti g 
(6) 


II 


-0.221 


0.086 


0.547 


-0.664 


-0.203 


0.104 


III 


-0.215 


0.083 


0.529 


-0.646 


-0.199 


0.099 



that the t\ 

,2 (3p 



and t\ u energy spectra are the analogue of 
ID, 1 S) and p 3 ( 4 5, 2 D, 2 P) terms in atomic 
This occurs because p x , Py and p z orbitals 
and Yj°) also belong to the t lu irrep of 1^ 



p 

physics 1 13 

fl32| . In Tabl e |IX| as well as in all other Tables of this sec- 
tion (X-XIV) the energy associated with the spherically 
symmetric multipole component (i.e. (2)^0 f° r electrons, 
( 2~ m )^o for holes and Ae for excitonic configurations) 
is put to zero. 

Similarly, one can obtain the energy levels of t\ g and 
Table |xj. Notice that the energy span of t\ u and t\ u 
configurations is almost the same, ~ 0.5 eV, and it is 
smaller than that of t\ g and t?„ states, ~ 0.75 eV. How- 



/ 3 



lg 



ever, the zero of energy in Tables |ix| and |x| is different 
for t™ u and £™ configurations. When one t\ u electron is 
promoted to a t\ g state, its energy is increased by ^1.153 
eV P, that is, 



Aei = e(h g ) - e(t lu ) » 1.153 eV. 



(5.5) 



This one-electron energy difference accounts for the in- 
teraction of the electron with the carbon nuclei and the 
"core" like <r— and it— electrons. Therefore, comparing 
the energy of {t\ g ) 2 with that of (ii„) 2 states one should 
add 2Aei to the (ii g ) 2 values. For the case of three 
electrons ((ii«) 3 and (ii g ) 3 ) we add 3Aei to the (ii g ) 3 
values. 

Next, in Tables XI and [Xll| we give results for the hole 
configurations (h^) m . We observe that the energy span 
of (h+) 2 and (/1+) 3 , ~ 1.2 eV, is almost the same. The 
magnitude is larger _than for_electronic t\ u and t" g con- 



figurations, Tables IX and [x|. Even a larger value of 
energy splitting, ~ 2.4 eV, was obtained for the case of 
four and five holes, (h^) 4 and (h+) 5 . Another important 
observation is that the number of states in a small en- 
ergy interval Ae ~ 0.03 eV near the ground state is 30, 
40, 25, 6 for m =2, 3, 4, 5, respectively. This suggests 
that the configuration of m = 3 holes is most susceptible 
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TABLE XI: Molecular terms of {hi) 2 and (h+f calculated 
with models II and III; in eV. ( 10 2 -m )£/o is zero of energy. 









(hi) 


3 






dee 1 


II 


III 






deg. 


ii 


III 


3 Tin 


(9) 


-0.248 


-0.242 




i lu 


(12) 


-0.707 


-0.687 


3 G„ 


(12) 


-0.223 


-0.216 




Gu 


(16) 


-0.682 


-0.662 




(9) 


-0.219 


-0.213 




4m 
J-2u 


(12) 


-0.678 


-0.658 


1 G g 


(4) 


-0.047 


-0.049 




2 rri 

-L2u 


(6) 


-0.430 


-0.423 


l H g 


(5) 


-0.040 


-0.042 




2 H 

u 


\ lu ) 


-0.414 


-0.407 


l H g 


(5) 


0.410 


0.402 




2 T 
J- lu 


(6) 


-0.401 


-0.394 


X A a 


(1) 


0.989 


0.958 




2 H U 


(10) 


-0.396 


-0.389 












G u 


(8) 


-0.077 


-0.075 












u + T'2u 


(12) 


-0.065 


-0.062 














(10) 


0.156 


0.155 












G u 


(8) 


0.160 


0.159 












2 H U 


(10) 


0.531 


0.514 



TABLE XII: Molecular terms of (hi) 4 and (hi) 5 calculated 
with models II and III; in eV. ( 10 ^~ m )?7o is zero of energy. 



(h+) 4 


(hi) 5 




deg. 


II 


III 




deg. 


II 


III 


b H g 


(25) 


-1.377 


-1.337 


6 A 

/lit 


(6) 


-2.294 


-2.228 


G g 


(12) 


-1.005 


-0.982 


H u 


(20) 


-1.716 


-1.672 


\H g 


(15) 


-1.004 


-0.981 


H u 


(10) 


-1.324 


-1.294 


G g 


(4) 


-0.819 


-0.805 


Gu 


(8) 


-1.323 


-1.293 


1 A 


(1) 


-0.814 


-0.799 


4 H U 


(20) 


-1.265 


-1.228 


3 T lg 


(9) 


-0.787 


-0.767 


Gu 


(16) 


-1.258 


-1.221 


Gg 


(12) 


-0.780 


-0.760 


2 A 
/lit 


(2) 


-1.106 


-1.079 


li g 


(9) 


-0.766 


-0.747 


Gu 


(8) 


-1.094 


-1.067 




(5) 


-0.491 


-0.482 


12u 


(12) 


-1.086 


-1.057 


T 2 g 


(3) 


-0.490 


-0.481 


*~Ju 


(16) 


-1.083 


-1.054 


l H g 


(5) 


-0.570 


-0.460 


1 lu 


(12) 


-1.057 


-1.028 


'Tig 


(3) 


-0.461 


-0.451 


2 H U 


(10) 


-0.816 


-0.800 


3 rri 

li g 


(9) 


-0.445 


-0.430 


2 Ti u 


(6) 


-0.803 


-0.787 


3 H g 


(15) 


-0.440 


-0.426 


2 H U 


(10) 


-0.792 


-0.775 


3 Ti g 


(9) 


-0.434 


-0.419 


2 T2u 


(6) 


-0.771 


-0.755 


3 H g 


(15) 


-0.434 


-0.419 


2 H u 


(10) 


-0.769 


-0.748 


G g 


(12) 


-0.160 


-0.157 


2 T2u 


(6) 


-0.769 


-0.747 


3 T2 9 


(9) 


-0.160 


-0.157 


2 H U 


(10) 


-0.747 


-0.726 


3 Ti g 


(9) 


-0.150 


-0.148 


2 T\u 


(6) 


-0.740 


-0.718 


1 A 


(1) 


-0.148 


-0.143 


2 Tvu 


(6) 


-0.487 


-0.476 


G g 


(4) 


-0.138 


-0.134 


Gu 


(8) 


-0.453 


-0.442 


l H g 


(5) 


0.020 


0.014 


2 T2u 


(6) 


-0.432 


-0.421 


1 G 9 


(4) 


0.041 


0.035 


Gu 


(8) 


-0.417 


-0.401 


\H g 


(5) 


0.087 


0.089 


2 A 

flu 


(2) 


-0.411 


-0.395 


G g 


(4) 


0.090 


0.091 


Gu 


(8) 


-0.231 


-0.223 


l H g 


(5) 


0.484 


0.471 


2 H U 


(10) 


-0.227 


-0.218 


l A a 


(1) 


1.062 


1.027 


2 H u 


(10) 


0.144 


0.135 



for Jahn- Teller distortions of the Ceo molecule and hence 
for hole-phonon coupling which causes superconductivity 

The results of calculations of excitonic confi gurations 
(/i+ij~ u ) and (h+ti) are quoted in Table XIII . The en- 
ergy span of (hit^g), ~ 1.5 eV, greatly exceeds that of 
(hiti u ), ~ 0.57 eV. A promotion of one electron to the 



TABLE XIII: Molecular terms of the excitonic configurations 
(hit^ u ) an d (h tt^ ) calculated with models II and III; in eV. 
Ae in Eq. ( 4.4 ) is zero of energy. 



(ft£*r w ) 






deg. 


II 


III 




deg. 


II 


III 


A G g 


(12) 


-0.034 


-0.034 


J-2u 


(9) 


-0.244 


-0.240 


3 Ti g 


(9) 


0.014 


0.013 


3 T lu 


(9) 


-0.118 


-0.115 


3 H g 


(15) 


0.014 


0.013 


3 H U 


(15) 


0.117 


0.114 


3 T 2g 


(9) 


0.078 


0.076 


G u 


(12) 


0.123 


0.122 


G g 


(4) 


0.120 


0.112 


1 H U 


(5) 


0.208 


0.200 


'Tig 


(3) 


0.187 


0.178 


G u 


(4) 


0.247 


0.238 


X T2 9 


(3) 


0.197 


0.189 


Tzu 


(3) 


0.513 


0.495 


"H g 


(5) 


0.548 


0.534 




(3) 


1.341 


1.321 



TABLE XIV: Molecular terms of (t lu t lg ) and (t lu ) 2 t lg cal- 
culated with models II and III; in eV. 2 G g and 2 ?2 S are "ac- 
cidentally" degenerate. (2)^0 is zero of energy. 





(t 


Lutlg) 






(tluf 


tig 






deg. 


II 


III 




deg. 


II 


III 


-flu 


(1) 


-0.969 


-0.958 




(20) 


-1.028 


-1.010 


3 H U 


(15) 


-0.474 


-0.467 


2 Ti g 


(6) 


-0.991 


-0.979 


3 7^ 

J lu 


(9) 


-0.390 


-0.384 


4 A 


(4) 


-0.901 


-0.886 


1 H U 


(5) 


0.473 


0.465 


Gg + l2 g 


(14) 


-0.398 


-0.395 


3 A 


(3) 


0.503 


0.506 


4 T lg 


(12) 


-0.376 


-0.361 




(3) 


0.697 


0.682 


2 H g 


(10) 


-0.223 


-0.224 










2 Ti a 


(6) 


0.004 


-0.008 










2 H g 


(10) 


0.449 


0.443 










2 A g 


(2) 


0.729 


0.714 










2 Ti g 


(6) 


1.052 


1.045 



t\ u shell increases the one-electron energy by the factor 

Ae 2 = e(t lu ) - e(h u ) » 2.69 eV. (5.6) 

The quantity Ae 2 is called electron affinity of Cgo and it 
was measured experimentally || [3?} . It accounts for the 
energy difference due to the interactions of the electron 
with the carbon nuclei and the "core" electrons. The 



value should be taken into account in Eq. (4.4) (i.e. Ae 



Ae2) when (h^t lu ) is compared with the ground state 
energy of the neutral molecule. In the case of (hit^ g ) 
one should use Ae = Ae 1 + Ae 2 = 3.85 eV in Eq. (O). 



The molecu lar t erms for (ii u ii g ) and t\ g are 



given in Table XIV . The energy span of the excited con- 
figurations is relatively large. It is approximately 1.6 eV 
for (t\ u tig) and 2 eV for (ii„) 2 ii g . Since both configu- 
rations imply the excitation of one t\ u electron to a ti g 
state, we should add Aei to the energies of the (ti u tig) 
and (ti u ) 2 tig molecular terms, when we compare them 
with those of the (ti u ) 2 and (ti u ) 3 configurations. The 
(ti u ) 2 and (tiutig) (as well as (ti u ) 3 and (ti u ) 2 ti g ) groups 
of terms are of different parity and thus there is no config- 
uration mixing between them. Therefore, although some 
two-electron molecular terms of the (ti u ) 2 and (ti u ti g ) 
configurations overlap, they do not interact with each 
other. The same holds for the (ti u ) 3 and (ti 
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e(t 1fl )-e(t 1u ), eV 



VI. MAGNETIC MOMENTS 

In this section we will calculate the magnetic moments 
of Cg^ for different orientations of the molecule. In a 
small external magnetic field H we ad d to a many body 
Coulomb interaction V{r, r"), Eq. (3.6), a magnetic term 



V, 



mag 



= -M z ■ H, 



(6.1) 



where M. z = X)fc=i M z (k) is a sum of one-electron (one- 
hole) terms with 



M(k) = fi B (L(k)+2S(k)). 



(6.2) 



FIG. 3: Three lowest levels of the coupled (t lu ) 2 + (t lg ) 2 
configurations as a function of Aei = t(ti g ) — e(ti«). The 
ground state is the 1 A g singlet for Aei < 0.58 eV, and the 
[i T lg triplet for Aei > 0.58 eV. 



figurations. 

However, there can be a hybridization between terms 
of the same symmetry of (ti u ) 2 and (tig) 2 configurations. 
The (tig) 2 configuration requires promotion of two elec- 
trons to the tig shell, with a subsequent energy increase of 
2Aei ~ 2.3 eV. Since the value is relatively large, the hy- 
bridization is expected to be weak. In order to study this 
issue we have carried out calculations where the mixing 
between the (ti u ) 2 and (ti g ) 2 configurations was allowed. 
In the calculation we have considered couplings between 
two 3 Tig levels, two 1 H g levels and two 1 A g levels at dif- 
ferent values of Aei. As before, we have employed the 
method described in Sec. III. We have found that the en- 
ergy spectrum separates in two groups. A group at lower 
energies originates from the former (ti u ) 2 levels, while 
the other group at higher energies has a large parentage 
of the (tig) 2 states. In Fig. 3 we plot the energies of 
the three lowest levels as a function of Aei. An inter- 
esting feature of Fig. 3 is the crossing of the 3 Xi g triplet 
with the 1 A g singlet at 0.58 eV with subsequent inver- 
sion of their positions. Thus, if Aei < 0.58 eV then the 
ground state is the 1 A g singlet, while for Aei > 0.58 eV 
the ground state is the T\ g triplet. This unusual be- 
havior explains why the l A g singlet was reported as the 
ground state of by Negri et al, Ref. |l| (QCFF/tt 
method). From our calculation (Fig. 3) it follows that 
the reason for this is a small energy difference between 



Here \ib is the Bohr magneton, k = 1 — n for elect ron s 
and k — 1 — m for holes. The magnetic moment (6.1) 



tig and ti u states. In Ref. [L6| Aei=0.64 eV, which is onl 



half of the experimental value of 1.153 eV for C 60 
Such low lying ti g states lead to an overestimation of the 
1 Ag [t\ u ] — 1 A g [t\ g ] configuration mixing and lowering 
of the bonding 1 A g term below 3 Ti g . The experimen- 
tal value Aei = 1.153 eV || [38| implies that the ground 
state is a triplet, as obtained by our calculations and in 
accordance with Hund's rules. 



belongs to the class of one-particle operators discussed 
in Sec. III. Explicitly, for the two-particle case we find 

(I\Vmag\J) = [<tl|M,|ji) 5(i 2 j 2 ) + (i 2 \M z \j 2 ) 5(idi) 
-{h\M z \j 2 ) Sfoji) - (i2\M z \j x ) 5(iih)] ■ H, (6.3) 

where (i\M z \j) stands for the one-particle matrix ele- 
ments. The generalization of the procedure for a many 
particle case is given in Sec. III. 

The one-particle matrix elements of spin momentum 
are given by the standard expressions |34| . They are 
independent of the orientation of the Cgo molecule. In 
order to calculate the orbital polarization we start with 
the Ceo molec ule in the orientation of Cohan |32|. By 
means of Eqs. ( |2.2a -c) and taking into account that 



L Z Y^ C = imY^ s , 



L Z Y™' S = -imY l m ' c , (6.4) 



we obtain for the orbital momentum of t\ u states 

(ti w> 2|£,|tit,,3}=2.5» > (6.5a) 
(t lu ,3\L z \t lu ,2) = -2.5i. (6.5b) 

The other matrix elements of L z are zero. (Here and 
below all values of magnetic moments are given in /j,b-) 
In the case of h u states we find 

(h u ,2\L z \h u ,3) = ii, (h u ,3\L z \h u ,2) = ~i, (6.6a) 



(h u ,4\L z \h u ,5) = i, (h u ,5\L z \h a ,4) = -i, 



(6.6b) 



and the rest is zero. 

There are two equivalent approaches to study the Cgo 
molecule in the magnetic field. The first is to use an 
active operator R(lo), which rotates the molecule as a 
three-dimensional object. In such case the magnetic field 
is always directed along the z-axis, while the position 
of the molecule is specified by three Euler angles w = 
(a, (3, 7). In the second case the position of the molecule 
is fixed but the direction of the magnetic field is changed. 
In the latter case one has to know the matrix components 
of the three projections of molecular orbital momentum. 
Below we have adopted the first approach which is more 



familiar to us from our previous study of rotator functions 
fu] , [33[ . The advantage is that we are working only with 
the z-component of orbital momentum. The details of 
the technique are given in Appendices A and B. 

Having calculated the matrix elements of V mag as a 
function of the molecular rotation u, we diagonalize the 
matrix H = Y, a ,b v (^a,r b ) + V mag (a). The degeneracies 
of molecular terms are lifted and the magnetic moment 
of each sublevel p is given by 



0.55 -, 



M(p) = ( P \M z \p), 



(6.7) 



where \p) is the corresponding eigenvector. 

For two electrons (or two t\ u holes) we obtain M. z = 
(±4.5, ±2.5, ±2, ±0.5, 0) for 3 T lg (ground state), while 
M z ( 1 H g ) = (±5, ±2.5, 0). In the 1 A g state the spin 
and the orbital momenta are quenched and M z = 0. 
For three electrons we have A4 Z ( 4 A U ) = (±3, ±1) (the 
ground term); M Z { 2 H U ) = (±6, ±4, ±3.5, ±1.5, ±1) 
and M z ( 2 T lu ) = (±3.5, ±1.5, ±1). We immediately 
conclude that the coupling scheme of orbital and spin 
momenta is different from the atomic case. In order 
to clarify this issue we have studied the polarization 
of spin and orbital momenta s epa rately. By excluding 
the spin momentum from Eq. (6.2) we have found that 
C z ( 3 T lg ) = (±2.5(3), 0(3)), and C,Q-H g ) = (±5, ±2.5, 
0) (numbers in parentheses stand for degeneracy). In the 
molecular term 4 A u the orbital momentum is quenched, 
C z = 0(4), while C Z ( 2 H U ) = (±5(2), ±2.5(2), 0(2)), 
and £ z ( 2 T lu ) = (±2.5(2), 0(2)). By excluding the or- 
bital momentum from Eq. (p^2) we have obtained the re- 
sults expected from the spin multiplicity index of molec- 
ular terms: 2S Z = (±2(3), 0(3)) for 3 T lg (spin triplet 
state) and 2S Z = 0(5) for 1 H g (spin singlet). For t\ u 
we find 2S Z { 4 A U ) = (±3, ±1), 2S Z ( 2 H U ) = ±1(5), and 
2S z { 2 T lu ) = ±1(3). 

In Sec. V we have already discussed the effect of mix- 
ing between the (ti u ) 2 and (ti g ) 2 configurations on the 
energy of the ground state of the Cgp anion. The hy- 
bridization affects also the magnetic moments of the 3 T\ g 
ground state, which are given by 




FIG. 4: 5 for the ground state triplet Ti s as a function of 
Aei, C sn . The corresponding magnetic moments are given 
by Eq. (ill). 



TABLE XV: The magnetic moments M and the orbital mo- 
menta C of (hu) 2 , in hb- 





M g 




3 Tl9 


0; ±(l,3,4,5)g 


0.5 


0,0,0, ±(1,1,1)5 


Gg 


±(l,l,3,3,5,5) ff 


0.5 


±(1,1,1,1,1,1)5 


3 T 2 g 


0,0,0, ±(l,l,l)ff 


2 


0(9) 


Gg 


±(1,1)5 


0.8333 


M 


l H g 


0; ±(l,2)s 


0.0668 


M 


l H g 


0; ±(1,2) 5 


0.2335 


M 


1 A B 












reason behind this. We consider the result as unexpected, 
because the magnetic moment of a 7r or 5 MO of diatomic 
molecules is anisotropic in respect to the direction of the 
magnetic field. From our previous study of C^~ in a cubic 
environment it also follows that the crystal field of C^~ 
exhibits strong dependence on its orientation [l4| . In 
case of the Ceo molecule, the orbital t\ u (or h u ) space of 
the icosahedral symmetry is greatly reduced in compari- 
son with the 11 dimensional I = 5 space of the rotation 



M = 0, ± 5 , ±2, ±(2± 5 ), ±(4 + .9). 



(6.8) 



The magnetic moments of the unhybridized 3 T\ g triplet 
of the pristine t\ u configuration correspond to g = 0.5 
(dashed line in Fig. 4). The evolution of the (7-factor as 
a function of Aei is given in Fig. 4. 

The results fo r the hole configurations (h+) m are 



quoted in Tables XV - XVIII . In the case of four or five 



holes the number of molecular terms is too big (27) and 
we give only magnetic moments for ten lowest states. In 
general, magneti c m oments are described by an expres- 
sion of the type (|6.8| ) although in some cases two distinct 
values of g are required. 

Interestingly, we have found that the calculated mag- 
netic moments are independent of the molecular orienta- 
tion. The conclusion holds for both t lM and h u shells and 
we think that there must be a hidden group-theoretical 



TABLE XVI: Magnetic moments M of (h+) 3 , in \i B . 2 Ti 
and 2 T2 U states are degenerate. 





M 




5 


J- lu 


±(1,2,3,5,6,7)5 




0.5 


G u 


±(1,1,3,3,5,5,7, 7)5 




0.5 


ilu 


±(1,1,1,3,3,3)5 




1.0 




±(1,1,1)5 




1.0 


2 H U 


±(5i,52,l,2 -52,2 - 


5i) 


0.6529, 0.8265 


Tin 


±(5,1,2-5) 




0.25 


2 H U 


±(51,52,1,2-52,2 + 


5i) 


0.2218, 0.3891 


G u 


±(5,5,2-5,2-5) 




0.8983 


2 rri | 2 /T"l 

-tlti + J2u 


±(1,1,1,1,1,1)5 




1 




±(5i,52, 1,2 - 52,2- 


5i) 


0.9578, 0.9789 


G u 


±(5,5,2-5,2-5) 




0.5650 


2 H u 


±(0,1,2,3,4)5 




0.5 
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TABLE XVII: Magnetic moments M for ten lowest molecular 
terms of (/i + ) 4 , in hb- 





M 


g 




0, ±(1,2, 2, 3, 4, 5, 6, 6, 7, 8, 9, 10) g 


0.5 




±(g,g,2-g,2-g,2 + g,2+g) 


0.310 


3 H g 


0,±(gi,g2,2 - fl2 ,2 - <?i,2,2 + fll ,2- 


- g 2 ) 0.046, 0.092 


Gg 


±(1,1)3 


0.011 


1 A 








3 T lg 


0, ±(5,2-5,2,2 + 5) 


0.052 


Gg 


±(g,g,2-g,2-g,2 + g,2 + g) 


0.023 


3 T 2 g 


0,0,0, ±(1,1,1)5 


2.0 




0, ±(l,2) fl 


0.146 




0,0,0 






hole). The matrix L^j then reads 



Ln - 





Va 

2 



2 ' 







n x i 







\n z i 

^2 n V 1 
2^x * 





~2 n V 1 



I by L 

n-r i 





n, i 



2 '"V 1 

-n z i 







.(6.12) 



The matrix has the same magnetic moments (eigenval- 
ues), which are 0, ±1/2, ±1 (in hb), for any direction of 
H. The reasoning given above is suggestive and we are 
looking for a full-scale group-theoretical solution to this 
problem. 



TABLE XVIII: Magnetic moments M for ten lowest molec- 
ular terms of (/ij) 5 , in /jb. A4( 4 H) stands for ±((71,32, 1,2 — 
g 2 ,2-g 1 ,2 + g 1 ,2+g 2 , 3, 4 - 52, 4 - gi). 





M 






3 


6 A 


±(1,3,5)5 
M{ A H) 






1.0 








0.533,0.766 


2 H U 


±(31,32,1,2-52,2- 


3i) 




0.873,0.936 


Gu 


±(3,3,2-3,2-3) 






0.366 


A H U 


M{ A H) 






0.866,0.933 


G u 


±(3,3,2-3,2-3,2 


f 3,2 + 3,4 


-3,4 


-3) 0.167 


2 A 


±3 






1.0 


G u 


±(3,3,2-3,2-3) 






0.978 


4m 
J- 2u 


±(1,1,1,3,3,3)3 






1.0 


G u 


±(3,5,2-5,2-3,2 


4-3,2 + 3,4 


-3,4 


-5) 0.5 



VII. ELECTRON OPTICAL TRANSITIONS 

In this section we consider only the electronic dipo- 
lar transitions and the corresponding optical lines for 
Cg,7. The picture is not complete because there exist 
electron- vibration interactions ( "Herzberg- Teller" cou- 
plings) which can alter the symmetry of the initial or 
the final state Q. Here we omit the electron- vibration 
couplings and limit ourselves to the electronic part of the 
problem. 

The optically active transitions are associated with a 
nonzero expectation value of the electron dipolar opera- 
tor P. Since the electric-dipole moment, 



(7.1) 



group 15*0(3). However, this is not accompanied by an 
anisotropic behavior of orbital momenta. 

In order to understand this issue we have considered 
a simplified case of one electron on the t\ u molecular 
level. Applying H in a direction n = H/H, where n x = 
sinficos^, n y — sin Q sin (f> and n z — cosfi, we find 



Vmag — HL — H Lfi. 



Here the matrix Lf L is given by 



Ln — n x L x + TiyLy + n z L z . 



(6.9) 



(6.10) 



By using Eqs. (2.2a-c) for three t\ u MOs, after calculating 
the matrix elements, we arrive at 



—Mqu v i Mqti x i 
MgUy i M(jn z i 
—Mqji x i —MqUz i 



(6.11) 



where Mq = 2.5 [Ib- The magnetic moments are obtained 
through the diagonalization of L^. We find that M = 
and M. = ±2.5 ^b for any direction of H. The same 
conclusion is obtained for the case of one /i u -electron (or 



is an odd quantity in respect to the inversion symme- 
try, it follows that P has no matrix components between 
states of the same parity. Therefore, all spectral lines 
due to electric-dipole radiation arise from transitions be- 
tween states of opposite parity (Laporte's rule) |l3| and 
the following schemes are relevant for the Cg 
molecular ions: 



,; and Cg 



(t lu ) 2 
(hu) 3 



(7.2a) 
(7.2b) 



These configurations have been considered already in pre- 
vious sections, and now we can proceed to calculations 
of optical transitions. 

In atoms there are several additional selection rules 
which greatly facilitate line assignments. These rules are 
not developed for the icosahedral symmetry and in the 
following we have to rely on numerical analysis. The 
total intensity for the line from level A to level B is given 

by H 

I{A,B)=N(a)hv-^— S(A,B), (7.3) 

where N(a) is the number of CgQ~ molecules in state a, 
v = [Ea — Es)/h is the frequency and a = hvjc is the 
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wave number. Finally, S(A, B) is the line strength which 
is found as 



(7.4) 



ah 



The line strength is a very convenient quantity and in 
t he fo llowing we calculate S(A, B) for the transitions 
(^2a|b). 

Th e matrix elements of the dipole operator for the case 
( 7.2a ) read 



{a\p\b) = j2m{i\p\j'){j'\b), 



(7.5) 



/./' 



where \I) and |J') are the basis states of (t\ u ) 2 and 
{ti u tig), respectively, while (a\I) and (b\ J') are the eigen- 
vector s co rresponding to levels a and b. The dipole mo- 
ment (7A_) is a one-electron operator, Sec. III. Its matrix 
elements are given by 

(I\P\J>) = S(i ul ,j u ) (i u2 \P\j g ) - (i u i\P\jg) 5{l u2 , ]u ). 

We recall that (I\ = (« u i,i M 2|, where i u \ and i u2 (i u i > 
i u2 ) are indices referring to six t\ u states, and |J') = 
\j u ,jg), Eq. ( |3 . l[ ) . From parity consideration it fol- 
lows that the nonzero matrix elements are of the type 
(i u \P\j g ). In order to calculate them, we first rewrite P 
in the following form: 



Pr 



P„ 



Pz = 



'4tt 1c , 

Y rYl 

•Y°(f). 



(7.6a) 
(7.6b) 
(7.6c) 



Here Y{ are real spherical harmonics |3l[ , and as before 
r = \r |, while f stands for the polar angles (f2,</>). Then 
we find that the one-electron matrix elements of P are 



(iu\Pk\j g ) = Vc lyT(k) (i u ,j g ), 



(7.7) 



where r = (l,c), (l,s) or for fc = .T,y and z, re- 
spectively. The quantities ci = \, T (i u , j g ) are given by 
Eq. (3.10) for I = 1. In fact, these coefficients have been 
also used for the calculation of the (t\ u t\ g ) and (ti u ) 2 t lg 
configurations in Sec. IV and V. Finally, the radial part 



of Eq. (7.7) reads 



V 



J drr 3 K t iu(r)n tlg (r 



(7.8) 



Since we have already computed the eigenvalues and 
eigenvectors of (ti u ) 2 and (ti u ti g ) in Sec. V, we now ca n 



calculate the line strengths using equations (7.4)- (7 
The results are quoted in Table XIX 



TABLE XIX: Energies E ab = Aei + e ab (in eV) and line 
strengths (in V) of the transitions [(ti u ) 2 ; a] — ► [(iiuii s ); b], 
a = 1 — 3, b — 1 — 6, calculated with the model III. Only 
transitions with S(a, b) ^ are given. 





(tiuf; 


a = 3 T lfl 




1 A 


(tiutig); b 


tab 


S 


tab S 


tab S 












3 H U 


-0.350 


0.482 






J lu 


-0.267 


0.289 






1 H U 






0.392 0.482 




3 A 


0.624 


0.386 






Ti u 






0.609 0.161 


0.323 0.129 



TABLE XX: Energies E ab = Aei + e ab (in eV) and line 
strengths (in V) of the transitions {(tiu) 3 ; a} — * {(ti u ti g ; &}, 
a = 1 — 3, b — 1 — 10, calculated with the model III. Only the 
transitions with S(a, b) are given. 





(tiuf; 


a = 4 A u 


2 H U 


2 


Tiu 




tab 


S 


tab 


S 


tab 


S 
















2 Ti g 






-0.913 


0.002 


-1.103 


S < 10~ 3 


Gg 














2 f ' i 2rp 

Ij-u + 1 2u 






-0.330 


0.450 






4 Tlg 


-0.010 


0.771 










2 Hg 






-0.158 


0.120 


-0.349 


0.205 


2 Tlg 






0.058 


0.072 


-0.132 


0.120 


2 Hg 






0.508 


0.683 


0.318 


0.277 


2 A 9 










0.589 


0.193 


2 Tlg 






0.780 


0.602 


0.920 


0.362 



Similarly, one can treat the optical transitions (7.2b) 

o_ ^ 1 

for Cg . Now we consider the matrix elements of P be- 
tween three-electron basis states (I (ti u ) 3 \ = (11,12,23! 
and I J [(ti u ) 2 t lg \) = \ji,j2,j g ), and obtain 



(W) = 

(ii\P\j g ) S(i 2 ,ji)S(i 3 ,j 2 ) - (h\P\j g ) 8(i2,32)S(ia,3i) 
+ (i2\P\j g ) S(ii,j2)S(i 3 ,ji) - (h\P\3g) S(ii,ji)S(i 3 ,j 2 ) 
+ (i3\P\jg) S(ii,ji)5(i 2 ,j 2 ) - (i 3 \P\j g ) S(i 1 ,j 2 )6(i 2 ,j 1 ). 

Here again, the one- elect ron matrix elements (i u \P\jg) 
are specified by Eq. (7.7). The resulting line strengths 



are quoted in Table XX 
It follows from Table 



XIXl that for Cg there are three 

3n 



lines from the ground state T\ g , 

E( 3 T lg -> 3 H U ) = Aei - 0.350 eV, (7.9a) 

E( 3 T lg -> 3 T lu ) = Aei - 0.267 eV, (7.9b) 

E( 3 T lg -> 3 A U ) = Aei + 0.624 eV. (7.9c) 



With Aei = 1-153 eV, Eq. (|Jh, we obtain E( 3 T lg -> 
3 H U ) = 0.803 eV, E( 3 T lg -T^T^) = 0.886 eV and 
E( 3 T lg -> 3 A U ) = 1.777 eV. These values should be com- 
pared with two dominant bands at 1.305 eV (950 nm) 
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and 1.494 eV (830 nm) observed by the near- infrared ex- 
periments in solutions ||. We conclude that first two 
transitions can be tentatively ascribed to the experimen- 
tal bands if Aei is taken to be larger, Aei ~ 1.7 eV. 
Here it i s worth to notice that in our approach Aei in 
Eq. (|5.5| ) remains a phenomenological quantity which is 
not immediately connected with the term splittings. In 
Ref. Jl6| the authors have obtained that the ground state 
of C 6 q is the 1 A g singlet. We have shown in Sec. V that 
this is possible if Aei is small, see Fig. 3. However, if 
l A g is the ground state, then there is only one optical 



transition { l A g -> l T Xu at Aei + 0.323 eV, Table |XIX|) 
and comparison with the experiment becomes even more 
problematic. We conclude that our calculations are basi- 
cally in agreement with the experiment for G 2 ^ , although 
also a third band is expected. The position of the third 
band however can change due to the effect of configura- 
tion mixing discussed in Sec. V. 

The situation is less clear for the Cgg" molecular ion. 
Both our calculations and those of Ref. [t6] predict that 
the ground state is the A A U level. Then from Table XX 



we find that the only optical transition allowed by the 
selection rules is A A U — > A Ti g . This is at variance with 
the experimental consensus for three dominant bands at 
1350, 960 and 770 nm ||. In principle, the theoretical 
line A A U — ► A T\ g can be split by the crystal field and 
Jahn- Teller distortions, but the magnitude of the split- 
ting (~ 0.3 eV) seems excessive. The other possibility is if 
some transitions become allowed through the "Herzberg- 
Teller" (electron- vibration) mechanism ||. Further ex- 
perimental and theoretical investigations are needed to 
clarify the issue. 



VIII. DISCUSSION AND CONCLUSIONS 

We have presented an efficient configuration interac- 
tion method for many electron (hole) molecular terms of 
the molecular ion. The method is based on the mul- 
tipole expansion of the Coulomb interactions and takes 
into account the molecular symmetry. Although there are 
some similarities with approaches used for treating many 
electron effects in atomic calculations |39| , the technique 
is novel and original. Crystal electric field effects [Q and 
the spin-orbit coupling can be easily incorporated in the 
calculation. The technique can be used for other electron 
systems. 

We have applied the method for studying molecular 

terms of electron and hole configurations of Cg (n = 
2 — 4, m — 2 — 5), and for excitonic configurations (h^t^ u ) 
and (h^t Xg ) of the neutral molecule. In most of the cases 
the ground state is found to obey Hund's rules. 

Our calculations of the molecular term structure (t\ u ) 2 
for Cj^7 differs from the previous result of Negri et al, 
Ref. [lq. They have reported that the ground term of CggT 
is the~^A g singlet, while we have found that for realistic 
parameters (Aei = 1-15 eV ||, it is the 3 T\ g triplet, 
in accordance with Hund's rules. We have shown (Fig. 3 



and Sec. V) that the reason for this is that in Ref. [H] the 
one-electron energy difference between t\ g and t\ u states 
is too small, Aei =0.64 eV. 

Our results for the (h^) 2 hole configuration indicate 
that there are three very close (~ 0.03 eV) low lying 
molecular triplets of 3 2\ g , 3 G g and 3 T2 g symmetry. The 
lowest molecular terms for (h^) 3 belong to the A Ti u , A G U 
and A T2 U symmetry. The number of states in a small 
energy interval Ae ~ 0.03 eV near the ground state is 30, 
40, 25, 6 for hole configurations [h^) m with m —2, 3, 4, 5, 
respectively. This suggests that the configuration of m = 
3 holes is most susceptible for Jahn- Teller distortions of 
the Ceo molecule and possibly for hole-phonon coupling 
which causes superconductivity ||. 

We have calculated the magnetic moments of the (ti u ) n 
and (h+) m configurations, Sec. VI. The coupling of spin 
and orbital momenta differs from the Lande g— factor 
scheme of atoms. The magnetic moments do not depend 
on the orientation of the molecule with respect to an ex- 
ternal magnetic field. The latter statement was demon- 
strated explicitly for the case of one t\ u electron and one 
h u hole. We consider this as a group-theoretical puzzle of 
the icosahedral symmetry. We have also found new "ac- 
cidental" degeneracy between the 2 G g and 2 T2 9 states 



of t he (U u) tig and (ti g ) t\ u configurations (Tables |VII 



and XIV.) 



Finally, we have studied optical absorption associated 
with electron dipolar transitions (ti u ) 2 — ► (ti u ti g ) and 

(hu) 3 - 

lines ( 3 T lg 



{t\u) 2 ti g - For Cg we have found that two 



J iJ u and J Tj 



is 



Ti„) can be tenta- 
tively ascribed to the two near-infrared dominant bands 
at 950 and 830 nm. However, in addition a third band 
( 3 Ti g — ► 3 A U ) is expected from the calculation. For Cgp 
with A A U as the ground state, we have found that only 
one electron dipolar transition, A A U — > A T\ gi is allowed. 
It seems that a better understanding of optical transi- 
tions requires a study of the Herzberg- Teller (electron- 
vibration) effect which is beyond the scope of the present 
work. We suggest to perform optical experiments for Cgg7 
and Cg^ + in the gas phase to obtain more precise and 
full data on the optical lines which can shed light on the 
problem of electronic intra- molecular correlations. 
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APPENDIX A: L z FOR A ROTATED MOLECULE 



If the molecule is rotated away, then the MOs of t\ u 



and h u symmetry are given by Eqs. (2.2a-c) and (2.1a 
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e) in the coordinate system (x',y',z') attached to the 
molecule. The rotated functions can be expanded in 
terms of real spherical harmonics (RSH) defined in the 
fixed set of axes (x,y,z). For example, a rotation R 
defined by three Euler angles uj = (a, f3, 7) transforms 
ipi{hu) to 



L R(lo)Y 5 ° + J l~R(u)Y^ c . 
50 V 25 v ^ 5 



(Al) 



Here R(oj)Y£ and R(uj)Y^ c defines the rotations of F 5 ° 
and Y^ ,c , respec tive ly. T he rotated functions can be 
found from Eqs. ([f35|), and (B9), quoted in the Appendix 
B. Analogously we proceed with the other angular func- 
tions of tin and h u symmetry. In general, 



(A2) 



Here t, t' = (m, c), (m, s) or 0, and the rotator func- 
tions (matrices) U^, r (uA are specified in Appendix B, 
Eqs. (JB5|), (B7), and (B9). Now we are ready to cal- 
culate L z ip k (ti u ) (k = 1 — 3) and L z ip' k (h u ) (k = 1 — 5). 
By means of Eq. (6.4) we obtain 

(k\LJk')l. 



m>0 



s)k 9(m,c)k' 9(m,c)k 9(m,s) v), (A3) 



where the functions g T .k depend on u>, 

9(m,c)k(uj) = a{n) U[ m c)Ti {uj) + a(T 2 )U l { „ 



9(m,s)k 

(w) = a(n) U\ m s)Ti (uj) + a(T 2 )U\ m s)T2 {uj). 

Here a(r) stands for the coefficients of expansion of MOs 
of t\ u and h u symmetry in terms of RSH, 

^ fe (fi) =a(r 1 )Y l T '(n)+a(T 2 )Y l T2 (n), (A5) 

see Eqs. ( 2.2a -c) and ( 2.1a -c). For example, for the first 
MO (k = 1) of ti u symmetry we have T\ = 0, a{rx] = 
and r 2 = (5,c), a(r 2 ) = ^/7/25, and etc. The 



indices k,k' in (A3) belong to the same molecular shell 
(ti u or h u ), otherwise (fc|L z |fc')| w = 0. From Eq. (A3) we 
also conclude that 



(k\L z \k)\ u 
(k\L z \k')\ 



0, 

= (k'\L z \k)\: 



(A6a) 
(A6b) 



The former condition is a consequence of working with 
real spherical harmonics, the latter ensures the hermitic- 
ity of L z . 



For Y', we have 



Y'^ = RH Y,- m = J2 Y l lDl n-M. (Bib) 



Here D l nm stands for the Wigner functions given by 



(B2) 



d l (f3)nm is a reduced matrix element which is a real quan- 
tity, and Cnm = ±1 depending on n, m (see Eqs. (2.1.6) 
and (2.1.5) of Refs. |l]). From Eq. @ and the proper- 
ties 



d{P)- n ,-m = (-l) n+m d(P) nm , (B3a) 
C- n ,- m = {-l) n+m C nm = {-l) n+m C; im , (B3b) 



we find that 



BU-mH = DIM* 



(B4) 



We then combine ( Bla ) with ( Bib ) and use Eq. (B4) for 
deriving the transformation law of real spherical harmon- 
ics. After some algebra we find 



r>ym,c _ y o rjl 



+ 



U(n,c);(m,c) + Y l ' ^(n,s) ; (m,c) ^ ) 



n>0 



(B5) 



where c^ ;(mc) = V2ReD l 0m , and 



(n,c);(m,c) 



I 

(n,s);(m,c) 



-Im(D 



nm 1 n—m 



(B6a) 
(B6b) 



In Eqs. (B5)-( |B6b| ) and below for clarity we drop the 



argument uj. Analogously, rotating K ' we obtain 



ttl[ — X l tV .( m;j ) 



n>0 



where U 0;(m,3 



E \ Y l l C U {n,c);(m,s) 

V2 ImD l 0m , and 



X l U (n,s)-(m,s) 



Tjl 

(n,c);(m,s) 



1 

(n,s);(m,s) 



Jm(££ m - D< 
= Re(D l nm - D l n 



), (B7) 



(B8a) 
(B8b) 



APPENDIX B: ROTATION OF REAL 
SPHERICAL HARMONICS 

An active rotation R is specified by its Euler angles 
uj = (a,/3,7) pl| . It transforms a complex spherical 
harmonic Y™ to Y\ , where 

Y>™=R(u)Yr = Y, Y i nDl nm(")- 



(Bla) 



Finally, the rotation of Y® yields 

ry° = y? ui <0 + {yr u l M , + rr ^ s);0 ) , 



?i>0 



(B9) 



where U l . = D l 00 , U[ n c) . Q = s/2 ReD l n0 and 
-V2ImD l n0 . 



16 



T.L. Makarova, B. Sundqvist, R. Hohne, P. Esquinazi, Y. 
Kopelevich, P. Scharff, V.A. Davydov, L.S. Kashevarova, 
A.V. Rakhmanina, Nature 413, 716 (2001). 
J.H. Schon, C. Kloc and B. Batlogg, Nature 408, 549 

(2000) . 

J.H. Schon, Ch. Kloc and B. Batlogg, Science 293, 2432 

(2001) . 

K. Tanigaki and K. Prassides, J. Mater. Chem., 5, 1515 
(1995). 

O. Gunnarsson, Rev. Mod. Phys. 69, 575 (1997). 
M.S. Dresselhaus, J. Dresselhaus, and P.C. Eklund, Sci- 
ence of Fullerenes and Carbon Nanotubes, Academic 
Press, 1995. 

L. Forro and L. Mihaly, Rep. Prog. Phys. 64, 649 (2001). 
H. Kuzmany, B. Burger, and J. Kiirti, in Optical and 
Electronic Properties of Fullerenes and Fullerene-Based 
Materials, edited by J. Shinar, Z.V. Vardeny and Z.H. 
Kafafl, Marcel Dekker, 2000. 

C. A. Reed and R.D. Bolskar, Chem. Rev. 100, 1075 

(2000) . 

M. Boyle, K. Hoffmann, CP. Schulz, I.V. Hertel, R.D. 
Levine, E.E.B. Campbell, Phys. Rev. Lett. 87, 273401 

(2001) . 

R.C. Haddon, L.E. Brus, K. Raghavachari, Chem. Phys. 
Lett., 125, 459 (1986); R.C. Haddon, Science, 261, 1545 
(1993). 

D. A. Bochvar, E.G. Gal'pern, Dokl. Akad. Nauk SSSR 
Chem. 209, N.3, 610 (1973). 

E. U. Condon and G.H. Shortley, The theory of atomic 
spectra, (University Press, Cambridge, 1967). 

A.V. Nikolaev K. Prassides, K.H. Michel, J. Chem. Phys. 
108, 4912 (1998). 

A.V. Nikolaev and K.H. Michel, Eur. Phys. J. B 9. 



619 (1999); A.V. Nikolaev and K.H. Michel, cond 



mat/0112147 . 

F. Negri, G. Orlandi, F. Zerbetto, J. Am. Chem. Soc. 
114, 2909 (1992). 

R. Saito, G. Dresselhaus, and M.S. Dresselhaus, Chem. 
Phys. Lett. 210, 159 (1993). 

A. Ceulemans and P.W. Fowler, J. Chem. Phys. 93, 1221 
(1990); P.W. Fowler and A. Ceulemans, Mol. Phys. 54, 
767 (1985). 

A. Auerbach, N. Manini, and E. Tosatti, Phys. Rev. B 
49, 12998 (1994); ibid. 49, 13008 (1994). 



[20] J.L. Dunn and C.A. Bates, Phys. Rev. B 52, 5996 (1995) 

[21] M.C.M. O'Brien, Phys. Rev. B 53, 3775 (1996). 

[22] CP. Moate, M.C.M. O'Brien, J.L. Dunn, C.A. Bates, 

Y.M. Liu, and V.Z. Polinger, Phys. Rev. Lett. 77, 4362 

(1996). 

[23] L.D. Landau and E.M. Lifshitz, Quantum Mechanics - 

Non-relativistic theory, v. 3, (Pergamon Press, Oxford, 

1965), Chap. 102. 
[24] LB. Bersuker and V.Z. Polinger, Vibronic Interactions in 

Molecules and Crystals (Springer, Berlin, 1989). 
[25] I.D. Hands, J.L. Dunn, and C.A. Bates, Phys. Rev. B 

63, 245414 (2001). 
[26] B.N. Plakhutin, Chem. Phys. Lett. 227, 396 (1994). 
[27] J.M. Oliva, Phys. Lett. A 234, 41 (1997). 
[28] E. Lo and B.R. Judd, Phys. Rev. Lett. 82, 3224 (1999). 
[29] B.R. Judd and E. Lo, J. Chem. Phys. Ill, 5706 (1999); 

B.R. Judd and E. Lo, J. Phys. B 32, 1073 (1999). 
[30] B.N. Plakhutin and R. Carbo-Dorca, Phys. Lett. A 267, 

370 (2000). 

[31] C.J. Bradley and A. P. Cracknell, The Mathematical The- 
ory of Symmetry in Solids, (Clarendon, Oxford, 1972). 

[32] N.V. Cohan, Proc. Camb. phil. Soc. math. phys. Sci. 54, 
28 (1958). 

[33] K.H. Michel, J. R.D. Copley, D.A. Neumann, Phys. Rev. 

Lett. 68, 2929 (1992); K.H. Michel, Z. Phys. B Cond. 

Matter 88, 71 (1992). 
[34] M.Tinkham, Group Theory and Quantum Mechanics, 

(McGraw-Hill, New York, 1964). 
[35] W.I.F. David, R.M. Ibberson, J.C Matthewman, K. 

Prassides, T.J.S. Dennis, J. P. Hare, H.W. Kroto, D.R.M. 

Walton, Nature, 353, 147 (1991). 
[36] E. Tosatti, N. Manini, O. Gunnarsson, Phys. Rev. B 54, 

17184 (1996). 

[37] S.H. Yang, CL. Pettiette, J. Conceicao, O. Cheshnovsky, 

R.E. Smally, Chem. Phys. Lett. 139, 5 (1987); D.L. 

Lichtenberger, KW. Nebesny, CD. Ray, D.R. Huffman, 

L.D. Lamb, Chem. Phys. Lett. 176, 203 (1991). 
[38] T. Kato, T. Kodama, T. Shida, T. Nakagawa, Y. Matsui, 

S. Suzuki, H. Shiromaru, K. Yamauchi, Y. Achiba, Chem. 

Phys. Lett. 180, 446 (1991); M.A. Greaney, S.M. Gorun, 

J. Phys. Chem. 95, 7142 (1991). 
[39] I. Lindgren, and J. Morrison, Atomic Many-Body Theory 

(Springer- Verlag, Berlin, 1982). 



